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Abstract 

We extend in several ways a recently proposed method to construct one-dimensional 
chaotic maps with exactly known natural invariant measure {HQ]- First, we as- 
sume that the given invariant measure depends on a continuous parameter and 
show how to construct maps with robust chaos — i.e., chaos that is not destroyed 
by arbitrarily small changes of the parameter — and prescribed invariant measure 
and constant Lyapunov exponent. Then, by relaxing one condition in the approach 
of Refs. Ill 0], we describe a method to construct robust chaos with prescribed 
constant invariant measure and varying Lyapunov exponent. Another extension of 
a condition in Refs. [QIB] provides a new method to get robust chaos with known 
varying Lyapunov exponent. In this third approach the invariant measure can be 
computed exactly in many particular cases. Finally we discuss how to use diffeo- 
morphisms to construct maps with robust chaos, any number of parameters and 
prescribed invariant measure and Lyapunov exponent. 

Key words: nonlinear dynamical system, deterministic chaos, robust chaos, 
natural invariant measure 
PACS: 05.45.Ac, 05.45.-a 



1. Introduction 

The inverse problem for chaotic one-dimensional maps has been recently pro- 
posed and solved in some cases by Sogo [Q], 0]. Starting from a given invariant 
measure, the method allows to construct maps with that measure and Lyapunov 
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exponents in the form A = mm, with m = 2, 3, . . . This method can be used to 
construct exact examples of chaotic one-dimensional maps depending on the dis- 
crete parameter m and having exactly known invariant measure and Lyapunov ex- 
ponent. The first goal of this work is to extend this method to families of maps that 
depend on a continuous parameter, display robust chaos — i.e., they have a chaotic 
attractor which is not destroyed by arbitrarily small changes of the parameter — 
and have known invariant measure and Lyapunov exponent. 

Piecewise smooth maps may show robust chaos and have been used to de- 
scribe robust chaos in circuits yfl. On the other hand, many families of smooth 
maps have fragile chaos. For instance, the logistic map, x n+ \ = Arx n (1 — x n ), 
is chaotic for r = 1, but the attractor is periodic for a set of values of the param- 
eter r that is dense in the interval < r < 1 [4]. If such a dynamical system 
describes a real device, it may be impossible to know in advance whether the be- 
havior of the system will in fact be chaotic or periodic for some parameter value, 
which is always known with finite precision. However, robust chaotic behavior 
is required in many applications, including encrypting messages [HI, @], random 
number generators \m and engineering applications [8|]. Robust chaos also arises 
in neural networks [ as well as in the study of brain and population dynam- 
ics |Qt, 121. Thus, simple and easy ways of constructing robust chaotic attractors 
with exactly known properties, as those explored in this work, can be of interest 
in different fields of science and engineering. 

Andrecut and Ali first found a smooth map [13] and later a method of gen- 
erating smooth maps whose evolution is chaotic for whole intervals of the 
parameter. Another case is discussed in Ref. Ill5ll and we have recently explored 
several new ways to construct smooth one-dimensional maps with robust chaos 

The purpose of this work is to present other ways of constructing maps with 
robust chaos. The advantage of these new methods is that the Lyapunov exponent 
is known exactly and the invariant measure is known always with the first, second 
and fourth methods and in many particular cases in the third approach. The first 
three methods will be extensions of Sogo's inverse problem |li|2j]. In Sect. [2] we 
will show how to construct families of chaotic maps with a constant Lyapunov 
exponent and a given invariant measure depending continuously on a parameter. 
In Sect. [3] we present a method to construct families of chaotic maps sharing the 
same prescribed invariant measure and having a known Lyapunov exponent de- 
pending on the parameter. In Sect. |4] we will explore a method to construct robust 
chaos with a known Lyapunov exponent varying continuously with the parameter 
and a natural measure that can be computed exactly in many particular cases. Fi- 
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nally, in Sect. [5] and [6] we will discuss several ways to construct robust chaos with 
prescribed Lyapunov exponent and natural invariant measure by using diffeomor- 
phisms. 

We will consider one-dimensional maps on a finite interval [a, b), which for 
commodity will be reduced to [0, 1] by means of a linear transformation. 



2. Robust chaos with prescribed invariant measure 

We will use an easy extension of Sogo's method [Q], 0] to construct families 
of dynamical systems, x n+ i = f r (x n ), that are chaotic for a full range of the 
parameter r and have exactly known invariant measure and Lyapunov exponent. 
Let us assume that the maps f r are m-to-1 so that each value x G [0, 1] has m 
preimages y k , such that f r (y k ) = x, for k = 1,2, ... ,m. the natural invariant 
measure dp r = p r (x) dx satisfies the Frobenius -Perron equation |4|]: 

Prix) = f Pr (y)5 (x - f r (y)) dy = jr ^Mj. (1) 

To find solutions of this equation we will simplify it by further assuming that all 
the terms in the sum make the same contribution, so that the substitution yi — > x 
gives the condition 

Pr (fr(x)) \fr(x)\ = mp r (x). (2) 

Eq. © provides a useful practical way to construct a family of maps f r with a 
prescribed invariant density p r : just choose, for a given integer value of m > 2, a 
family of solutions of the differential equation (O that map the phase space [0, 1] 
onto itself and are m-to-1. This will be a family of maps with robust chaos. By 
extending the argument in Ref. [2] one can show that all the maps in the family 
will have the same Lyapunov exponent: 

X r = p r (x) In \fr(x)\ dx 
Jo 

— In m + / p r (x) In p r (x) dx 
Jo 

— / p r (x) lnp r (f r (x)) dx — mm. (3) 
Jo 
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The two last integrals cancels each other because of (OQ): 

Pr(x) In p r (x) d,X 

dyPr{y) dx 5 {x - f r {y)) In p r (x) 



p r (y) \np r (f r (y)) dy. (4) 

Jo 

2.1. Some solutions with m=2 

To give some examples, let us further assume m = 2 and that the invariant 
densities p r £ C°(0, 1) are strictly positive, p r (x) > for < x < 1, so that 

/x r (x) = / p r (y)dy (5) 



monotonously increases from p r (0) = to p r (l) = 1 and has a unique inverse 
p^ 1 {x) in [0, 1]. If we choose the boundary conditions / r (0) = f r (l) = 0, the 
solution of © is 

/ r (x)=/i T T 1 (l-|l-2/i r .(x)|). (6) 
Each map will increase from / r (0) = to the maximum 

.-i d 



/ r (a r ) = l, a r = p r I -J (7) 

and then decrease to f r .(1) = 0. The Lyapunov exponent will be A r = In 2. 

2.1.1. A family of piecewise smooth chaotic maps 
If one chooses the invariant densities 

p r {x) = rx r ~\ (r > 0), (8) 

the solution © will be 

/ r (x) = (l-|l-2x1) 1/r , (r>0). (9) 

Three members of this family of chaotic maps are displayed in Fig. [TJ The full 
tent map T(x) = 1 — |1 — 2x\ is recovered with r = 1. 
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Figure 1: Three members of the family of maps ©. 

2.1.2. A family of smooth chaotic maps 

The solutions © are not differentiable at x = a r . In fact, since 

l™ J'rix) = T 2 -^-, (10) 
x->c$ Pr{l) 

the condition for f r to be continuous at the maximum is lim^,^! p r (x) = oo. 

We can use this condition to find families of chaotic maps that are smooth 
along the full interval [0, 1]. To provide an example, let us consider the family of 
densities given for every r > by 

2 r ~ 1 r arcsin r_1 Jx~ 

» w= «v«Q-«) " <H) 

Then the solution © is 



sin 2 (2 1/V arcsin y/xj , < x < a r ; 

fr{x) = { . 2 f l/r{^ ■ r MlM . . , (12) 

1 sin I 2 i/r — arcsin y/x) J , a r < x < 1, 



with a r = sin 2 (2-( 1+1 / r )7r) . 

Some members of the family (fT2l) are displayed in Fig. [2l where one rec- 
ognizes the full logistic map fi(x) = 4x(l — x). We have f r 6 C^O, 1] for 
< r < 1, but ~~ * — 00 as x — ► 1 for r > 1. 
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Figure 2: Three members of the family of maps (fT2l . 

The boundary conditions for f r (x) can be chosen in many other ways. For 
instance, with / r (0) = / r (l) = 1 the maps will have a single minimum, instead 
of a maximum, and the substitute for © is 

f r (x)=^ l (\l-2 l i r (x)\). (13) 

The single minimum is f r (a r ) = with a r = /i" 1 (1/2). 
With the densities ® one gets 

f r (x) = \l-2x r \ 1/r , (r>0), (14) 

which has a smooth minimum only for < r < 1, since now the condition is 

]im x ^ p r (x) = oo. 

It is also easy to extend the results of this section for values of m ^ 2. But let 
us present another way to construct robust chaos. 

3. Robust chaos with constant invariant measure 

Following Refs. |Q1 0], to get Eq. §2$ we assumed that all terms in sum (OQ) 
had the same value. It is obvious that this simplifying condition can be relaxed 
in many ways. Let us explore a simple one. One could construct a family of 2- 
to-1 maps, with the same prescribed invariant density p(x) but with a generalized 
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relation between the terms in sum CQ), 



P{VI) -r-£Hr, (r>0), (15) 



\f'r{yi)\ \fr(V2 

for instance. With this assumption, one would replace the constant multiplicity m 
of Eq. ©by 

t \ J 1 ?' 0<s<a r ; 
m r (x) = < r (16) 
I 1 + r, a r < x < 1; 

where x = a r is the value where the two preimages coincide: yi (a r ) = y 2 (a r ). 
The family f r would satisfy the following equation: 

p(f r (x)) \ffa)\=m r (x)p{x), (r>0). (17) 

For example, for any p G C°(0, 1) with p(s) > and = p(y) cfa/), the 
solution of O with / r (0) = / r (l) = is 

L-'^Kx)), 0<x<a r ; 

with 

a r = . (19) 

By using (flTT) and (fl9l) the varying Lyapunov exponent is readily computed: 

X r = [ p(x)ln\f'(x)\ dx = ln(l + r) — lnr. (20) 

Jo l + r 

This value starts from lim r ^ K = 0, increases until its maximum Ai = In 2 and 
then decreases towards Hindoo A r = 0. The condition for the maximum at x = a r 
to be smooth is again lim^i p(x) = oo. 

3.1. Another family oj piecewise- smooth chaotic maps 

For example, if one selects for all r > the same constant invariant density 
p r (x) = p(x) = 1 the solution of Eq. dT8b is the family 

„ I x/r, < x < a r : 
f r {x) = {l + r)\ ' - - r : r>0, (21) 

I 1 — x, a r < x < 1; 

with a r — r/(l + r) and the Lyapunov exponent of Eq. (|20l) . /i is the full tent 

mapT(x) = 1 — 1 1 — 2x\. 
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3.2. A second family of smooth chaotic maps 

If the starting point is the invariant density of the full logistic map, p r (x) 

p(x) = [ir 2 x(l — x)] , we get, for every r > 0, 



sin ((1 + r) arccos y/x) , a r < x < 1; 



a r = cos2 TT7T~ — V ( 23 ^ 
2(1 + r) 



with 



and the Lyapunov exponent (|20l) . The full logistic map is recovered with r = 1: 

= 4x(l — x). 

It is straightforward to extend the method in this section for m = 3, 4, . . . 

4. Robust chaos with varying Lyapunov exponent 

In Eq. © it was assumed that m is an integer larger than one. But let us 
substitute for it a real number r > 1: 

ptf r {x))\f r (x)\=rp{x), (r>l). (24) 

Notice that we are using again the same invariant density p(x) for all values of the 
continuous parameter r, which now replaces the constant multiplicity m. If we 
get a solution of Eq. (|24l) for a given p(x), the latter will not satisfy the Frobenius- 
Perron equation for non-integer values of r, but this does not prevent the solution 
f r from being chaotic for all r > 1, as we will show in the following. In fact, by 
using an argument very similar to the one leading to ©, one can show that the 
Lyapunov exponent of the family satisfying (|24l) will be A r = lnr, for all r > 1. 

The solution of Eq. (gS) for / r (0) = 0, /;(0) > and 1 < r < 2 can be 
written as 

f ( x j = MX))> < x < a r = p^il/r); ^ 
1 p~ x (2 — rp(x)) , a r < x < 1, 

with p(x) = dy. It is also easy to write down the solution for other initial 

conditions or other ranges of the parameter r. 
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Figure 3: Bifurcation diagram of the map d26b . 



4.1. A third family of smooth chaotic maps 

If we choose again the invariant density of the full logistic map, p(x) = 
[ir 2 x(l — x)] -1 / 2 , a solution of Eq. (|24|) for all r > 1 is 

f r (x) = sin 2 (r arcsin yfx) . (26) 

For r = 2, 3, ... the maps f~ (x) reduce to polynomials and we recover the 'Cheby- 
shev hierarchy' of Refs. Ill, 120 : 

fn(x) = \- ^f- T n (2x - 1), (n = 2, 3, . . .), (27) 

where T n are the Chebyshev polynomials of the first kind. All these maps are 
chaotic, with Lyapunov exponent A n = In n and the same invariant invariant den- 
sity as the full logistic map f2(x) = 4x(l — x). The question is what happens for 
non-integer values of r? 

Their Lyapunov exponent is A r = In r and we can see in the bifurcation dia- 
gram of Fig. |3]that the attractor only fills the phase space [0,1] for r > 2. 

But, what about the invariant measure? The point is that now all points x do 
not have the same number m of preimages y k , such that f r (y k ) = x, as assumed 
in Refs. [Ql 0] and Sect. [2] and [51 In consequence, the starting invariant density 
will not satisfy the Frobenius-Perron equation (OQ) for non-integer values of r. 
Obviously, the actual invariant measure will satisfy that equation with a value of 
m depending on the point x. It it clear from Fig. @] which display in solid line 
the graphs of two members of family (|26l , that the value m(x) will change in this 
example at point x = a = f r (l) and the same will happen at all its images, so 
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Figure 4: The maps ([26) and OH for r = 6. 1 + \/2. 



that we can expect the invariant measure to be discontinuous at every point in the 
orbit of / r (l): 

O r ={f k r {l): k = 1,2,...}. (28) 

The task of finding the invariant measure can thus be very difficult unless the 
orbit O r is simple enough (a cycle), as happens in the 'Chebyshev hierarchy' 
(where r = 2, 3, . . . and 2n = {0, 0, . . .} and 2 n-i = {1,1,.. •}, so that the 
invariant density is continuous for < x < 1), but also for many non-integer 
values of r as we shall see now. 

For instance, in the left graph of Fig. |4] the parameter equals the golden ratio, 
r = = (1 + \/5)/2, and then the orbit of /^(l) is a 3-cycle: 

0$ = |a = sin 2 1 — a, 1, a, 1 — a, 1, . . . j> . (29) 

Since one can multiply a solution of Eqs. (Q]) and (l24l) with any constant, one may 
suspect that the actual invariant measure of the maps (|26|) is that of the full logistic 
map multiplied by a function which is constant except at the points lying on the 
orbit of f ^(l). In fact, we can find the exact expression of the invariant density 
for r = by using the normalization condition p^x) dx = 1 and trying in the 
Frobenius-Perron equation (OQ) a density in the form 

{a, < x < a; 
b, a < x < 1 - a; (30) 
c, 1 -a < x < 1; 

with constant a, b and c and the following multiplicity: 

1, 0<x<a; 
m{x) = { n (31) 
, a < x < 1. 
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One finally finds 

a = 0, b = — = . (32) 



The result can be easily checked by means of numerical simulations [17]. 
Even simpler is the orbit of / r (l) for r = 1 + y2: 

1+v /2 = ja = cos 2 -^=, 1, a, 1,...J. (33) 

From this 2-cycle and the Frobenius-Perron equation (0Q) it is easy to find the 
invariant density: 

2 + V2 {y/2, 0<x<a; 
Pl+V ^ 47iy/x(l-x) [1, a<x<l. 

Other cases with piecewise continuous invariant measure can be found in a 
similar way. 

It should be noticed that maps (1261 ) are a nice example of solvable robust chaos, 
since the solution can be explicitly written in terms of the initial condition xq as 

x n = sin 2 (r n arcsin y/xoj . 

5. Robust chaos through diffeomorphisms 

Another way to construct examples of robust chaos takes advantage of the 
invariance of Lyapunov exponents through diffeomorphisms on the interval [0, 1]. 
We are going to consider three different approaches. 

5.1. A single diffeomorphism 

If we know a family of chaotic maps f r (x) with invariant densities p r (x) and 
choose a diffeomorphism ip(x), the topologically conjugated family f r = ^p o fr° 
1 will have the same Lyapunov exponent and the following invariant density 

wvp 0*0)1 

For example if f r is the family (1211) and the diffeomorphism 

(p(x) = sin 2 (36) 
the topologically conjugated family is (|22l) . 
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On the other hand, if f r is the family (|26l) and the diffeomorphism 

2 

tpix) — — arcsin yfx, (37) 

7T 

the topologically conjugated family is 

f r (x) = 1 - |1 - rx mod 2| . (38) 

Maps (|38l) are piecewise linear and two particular cases are depicted in dashed line 
in Fig. HI Moreover, one recovers the full tent map for r = 2. It is obvious that 
since the slope of the graph of f r (x) is nearly everywhere ±r the maps are chaotic 
for all r > 1 and that the Lyapunov exponent is A r = lnr. This in turn proves 
again that the topologically conjugated family (fT2l) is chaotic and has the same 
Lyapunov coefficient. The invariant density is p n (x) = 1 for r = n = 2, 3, . . . 
and can be computed for other values of r by using (1351) or the method discussed 
in the last section. For example, the measure corresponding to (|34|) is piecewise 
constant: 

2 + V2 jV2, 0<x<V2-l; 
= — ^ — i ^ 1 ( 39 ) 



[1, V2-l<x<l. 



5.2. A family of diffeomorphisms 

A second possibility is to choose a single chaotic map f(x) with known in- 
variant density p(x) and Lyapunov exponent A and a family of diffeomorphisms 
depending on a continuous parameter: (f r {x). Then the maps f r = (p r o / o tp~ l 
will have the same Lyapunov exponent and the following natural density: 



«*) = WO) 

For instance, if the starting map is the full tent map, f(x) = T(x) = 1 — |1 — 
2x\, with p(x) = 1 and A = In 2 and the diffeomorphism (p r (x) = x 1//r , with 
r > 0, the family f r is precisely © with the same Lyapunov exponent and the 
invariant density ©. (Notice that in this example / = f\.) 

On the other hand, given the full tent map and the diffeomorphisms (f r (x) = 
sin 1 / 1 * with r > 0, the topologically conjugated smooth maps 

f r {x) = sin 1/r (2 arcsin x r ) = x (4 - 4x 2r ) 1/2r (41) 
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have the same Lyapunov exponent and the natural densities 



Pr(x) = )^= . (42) 



X 



This map, which is also obtained from the full logistic map f(x) = 4r(l — x) by 
means of the diffeomorphism (p r (x) = x x l 2r ', is a simple example of exactly solv- 
able robust chaos, in which everything is known, including the general solution 
x n = sin 1//r (2™ arcsin Xq). We recover the full logistic map with r = 1/2. 

5.3. Prescribed invariant density 

Finally, one can choose a chaotic map f(x) with known invariant density p(x) 
and Lyapunov exponent A and a prescribed family of natural densities p r (x) de- 
pending on the parameter r. Then each family of diffeomorphisms (p r satisfying 
the differential equation 

p r (iPr(x)) \<f' r (x) \ = p{x) (43) 

will provide a family of maps f r = <p r o f o ip~ l with the same Lyapunov exponent 
and the desired natural measure. 

If we want (p r (x) to be increasing, the solution of (|43l) is 



ip r (x) = p r (p(x)) , (44) 



with 



p r (x)= / p r (y)dy, p(x) = / p(y)dy. (45) 
Jo Jo 

For example if f(x) is the full tent map and p r (x) = rx r ~ l , the solution (1441) 
us (f r (x) = x x l T for r > 0, and we recover once more the family ©, which in 
turn shows again that the latter is chaotic for all positive real values of r and that 
its invariant density is ([8]). 

If we choose y? r (x) to be decreasing, the solution of (1431) is ipjx) = fx" 1 (1 — p(a 
and for p r (x) = rx^ 1 and r > we have tp r (x) = (1 — x) 1//r , which transforms 
the full tent map into the family (|T4|) . Since families (TT4)) and © are topologically 
conjugated to the full tent map, they are also conjugated to each other. In fact, the 
diffeomorphism conjugating them is (p r {x) = (p~ x (x) = (1 — x r ) 1 ^ r . This a nice 
example of two topologically conjugated maps sharing not only the Lyapunov 
exponent but also the natural invariant measure. 

This simple example shows that the maps found in previous sections are not 
necessarily the only solution to the corresponding problem. It was pointed out in 
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Refs. [Q], 0] that in general the inverse problem does not have a unique solution 
and several examples were described with different maps sharing a given invariant 
density while having different Lyapunov exponents in the form \ m = In m. We 
see here that it is also possible to have different maps with the same Lyapunov 
exponent and invariant measure. 

An obvious variant of the method in this section is a starting family of densities 
p r (x) depending on the parameter r to construct a family of chaotic maps with a 
constant density p(x) by solving 

p(tp r (x))\tff r (x)\ = p r (x) . (46) 



6. The inverse problem 

Obviously, one can combine the methods in the previous sections to get fami- 
lies of maps depending on several parameters which display robust chaos. Instead, 
we will discuss a more general and easier method to accomplish the same. 

All the maps considered before can be written in the form / = p^ 1 o g o p, 
for some appropriate g and p. In fact, for the particular case of smooth unimodal 
families with robust chaos it has been proven that all maps within the family are 
topologically conjugate [18]. This suggests a general way to construct a map 
/ (smooth or not, unimodal or not) with prescribed Lyapunov exponent A and 
natural invariant density p(x), such that p(x) = J* p(y) dy is a diffeomorphism, 
even in the case in which /, A and p depend on one or several parameters, which 
will not be written explicitly. 

The starting point is a set of m values < a k < 1 such that 

m 

b m = J2a t = l. (47) 

k=i 

If bo = and b k = + a&, for k = 1, 2 . . . , m, the piecewise linear map 

( X - htl h <T <b ■ 

9(x) = \ h ^~ X x (48) 
— , hi-i <x<b 2i 

(with i = 1, 2, . . .) has the invariant density po(x) = 1, because of (|47l) . and the 
Lyapunov exponent 

m 

< Ao = — a k In a k < In m. (49) 

k=l 
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The maximum Lyapunov exponent A = In m is obtained with a± = a 2 = • • • = 

a m = 1/m. 

Now, if one chooses a value of m large enough (so that A < In m) and a set 
of values a& (depending on the parameters of the desired map) such that A = Ao, 
then the map / = /i -1 o g o fj, will have the prescribed Lyapunov exponent A and 
natural invariant density p(x) = p'(x). The derivative /' will be continuous at 
the maxima (minima) if lim^x = oo (linx^o = oo). The map / will not be 
unique, since many others can be constructed in the same way with higher values 
of m and, even for the same m, there will be in general many (infinite) ways of 
choosing the values a^. For instance, this method reduces to (fT3l) for m = 2, 
o-i = a 2 = 1/2 and g(x) — T(x) = 1 — |1 — 2x\. On the other hand, solution ( TT81) 
is recovered with m = 2, a 2 = 1 — a x = 1/(1 + r) and g given by (I2TT) . 

It is also easy to change slightly this method to construct maps with known 
Lyapunov exponent and piecewise natural density. Let us consider a single case. 
Instead of (|48T) the starting point will be 



for m = 3, 5, . . . Since the orbit of g(l) is O = {a ly 1, ai, . . .} there will be a 
single discontinuity point in the natural invariant density, which in fact is 




(50) 




< x < a\\ 
a\ < x < 1, 



(51) 



while the Lyapunov exponent is 



in 




< A 



k=l 



< arcsinhn, n = (m — l)/2. 



(52) 



o — 



1 + a 



"in 



The maximum value is reached with a\ = ■ ■ ■ = a m _i = ^fa^ = y/i + n 2 — n. 
Then for each choice of the values a/-, the map / = pT l o g o jj, will have the 
Lyapunov exponent A = A and the natural invariant density 




< x < /i -1 (d) ; 

(Ol) < X < 1. 



(53) 
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For instance, taking m = 3 and a x = a 2 = \J2 — 1, one recovers the map (1381) 
for r = \/2 + 1, which can be extended to a piecewise linear family with known 
Lyapunov exponent and invariant measure by letting a 2 run from to 2 — y/2. In 
turn, applying to that family the diffeomorphism fj,(x) = - arcsin y/x one gets a 
family of smooth maps with known Lyapunov exponent and piecewise continuous 
natural invariant density, which coincides with (|26l) when a 2 = l/r = \[2 — 1 . 

It is easy to construct other examples starting from a piecewise linear map 
with simple orbits of g(l) (or g(0)). 

7. Final comments 

We have presented some new easy methods to construct families of one-dimensional 
maps with robust chaos and exactly known Lyapunov exponent and natural invari- 
ant measure. As far as we know, this is the first time general methods to do that 
are discussed. 

It should be stressed that the explicit examples with robust chaos discussed 
above have been selected for simplicity, but many other can be easily constructed 
with the ideas discussed in this work. 
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